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MECHANICS  AND  THERMODYNAMICS  OF  A  MIXTURE  OF  A 
GRANULAR  MATERIAL  WITH  A  FLUID 

S.  L.  Passman 


Introduction 

In  a  previous  work  [  1]  balance  equations  for  a  mixture  *>f  an 
arbitrary  finite  number  of  granular  materials  has  been  given;  however, 
no  constitutive  theory  has  been  developed.  In  this  work  I  consider  the 
special  case  of  two  materials ,  one  a  granular  material  as  defined  by 
Goodman  and  Cowin  f  2] ,  the  other  a  viscous  fluid.  A  constitutive  theory 
is  postulated,  and  restrictions  due  to  the  entropy  inequality  are  explored, 
following  closely  the  analysis  of  Muller  [  3] . 

I  use  without  further  comment  the  equations  and  notations  of  f  1] „ 
Furthermore,  sections  and  equations  herein  are  numbered  as  if  this  work 
were  a  continuation  of  [  l'j . 


5.  Further  Analysis,  of  the  Entropy  inequality 

Recall  the  equaMo.is  for  balance  of  energy  and  entropy  for  the 
mixture : 


pc  =  tr^graa  k)  +  h  *  grad  v  +  pfc*  2 


+  pgv  +  div  3  +  ps  , 


(5.1) 
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"1  ■ 


. . . 


jiaasiiiaa 


div  £  =  div  £  +  3  *  grad  3  +  3  div 

-  nV  grad  3  -  3Zu  •  div  T  -  32  tr(TTgrad  u)  (5,13) 

a  0  a  0  a  a 

-  2 h(v-O) .  grad ’3-'3Z^v -v)div  h  -  $h  •  grad(w-C')  . 

0  a  0  0  a  a 


#1391 


-3- 


Also,  the  Helmholtz  free  energy  of  the  mixture  is 


=  < 


15.14} 


Eliminating  b  and  f  from  (5.9),  and  using  (5.13)  and  (5.14) 
o  a 


gives,  by  use  of  (2.  25)^  , 


pn  1  -  p4  >  div  J  -  3[  2px  ‘  u  +  2p(k  v  +  icv){v  -  &)] 

ca  a  e  a  a  oo  o 

T 

+  tS  '  H  -  2h(w-  v) }  .  grad  5 

C  Q  Q  0 

t  + x  41  -t*  v  v 

+  tfp2(m  -  cx)  ♦  u  +  -3p2(v -  ckv)(v-  v) 
ooo  a  aoooo 


(5.15) 


18 1  tr(T  grad  x)  +  2  tr(T  grad  u)  +  h  *  grad  v  +  £h  •  grad(v-v  jJ 

q  o  ~  *  ii 

-  y'.pkv1  -2pgv)  . 

0  30 


Define 


G  =  grad  x  , 


(5.16) 


the  velocity  gradient,  and 


G  =  grad  x 

A/  M 

a  a 


(5.1?) 


the  peculiar  velocity  gracient  ci  the  a  -th  constituent.  The  symmetrl/and 
skew  parts  of  these  am 


(5.18) 


+  ST),  W  -  ~  (Q  -  GX), 

the  stretching  and  spin,  and 

5  =  ?(G+GT),  W=^(G-GT),  (5.19) 

a^’aa  a£oo 

the  peculiar  stretching  and  peculiar  spin  of  the  o-th  constituent. 

By  (2.11) 


grad  u  =  G-  G. 

*  •»  Pv  M  1 

a  a 


(5.2C) 


and  by  (1.6)  and  (1.10) 


so  that 


and 


pQ  =  ZpG  Uu®  grad  p  , 

a  a  o  q 


(5.21) 


PD  =  Sp  D  +  -  2(u®  grad  p  +  grad  p  &  u )  , 


a  a 


a  c 


(5.22) 


pW  =  2p  W  +  ~  2(u  ®  grad  p  -  grad  p®  u  )  . 


o  a 


a  a 


(5.22) 


By  (2.19)  and  (2.23)  ST  is  symmetric,  and  by  (2.25).  T  is 

^  Q  * 

symmetric.  Then  by  (5.18)  -  (5.21), 
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(5.23) 


tr(T  grad  x)  +  2  tr(TA  grad  u) 

a  a 

=  tr  TD  -  tr'Z  T)D  +  2)  tr  TTG 
a  a  a 

P  1  T 

=  tr(T  -  2)T)(E— •  D  +  ~  u  ®  grad  p)  +  2  tr  T  G 
a  o  P  a  Q  a  c 

The  gradient  of  (1.11)  is 

N  X 

(grad  p.k)v  +  pk  grad  v  =  r(grad  pk}v+2)pk  grad  v  . 

aa  g  oo  a 

However, by  (l.U), 


grad  pk  =  X  grad  p  k  . 

a  a 

Substituting  (5.  25)  into  (5.24).  I  obtain  the  result 


pk  grad  v  =  z  pk  grad  v  +  2[grad(pk)J{v  -  v)  . 

aa  a  ca  a 

Then 


h  •  grad  v  +  r  h  •  grad(v  -  v) 
a  o 

=  (h  -  2h)  «  •:}.-€}  4  v  +  Zh  •  grad  v 
a  oo 

=  (fc-  Sh)  *  '4; '  c.pk  grad  v  +  E  (v-v)grad(pk)] 
o  oo  a  o 

iih*  grad  '  . 
o  * 
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(5.24) 


(5.25) 


(5.26) 


(5.27) 
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Substituting  (5.23)  and  i5.2*  ?  into  (5. 15), yields 


^  A  ^  >  s  s 

pti  -  -  p3i  >  div  £  -  ’3[  I  px  •  u  +  E(pkv  +  pk^ )(v> - v ) ] 
^  aa  a  ggg  aaa  a 


+  (g  -  ETg-  Zh(’-'-v)]  •  grad  $ 
a  a  a  o 

"I'N  -f  ^  ^  9 

+  dp  Z(m-cx)  •  u  +  -3p  Z(v -ckv)(i'-v') 
aa  a  a  aca  a 

-  -  trfT  -  ETKEpD  +  Zu  ®  9rad  p)  -  3  Str  TT-G 

^  a  a  a  o  o  aa 

-  —  (h  -Eh)*  [  X  pk  grad  v  +  E(v  -v)(grad  pk)l 

P"^  n  ad  a  n  no 


(5.23) 


aa 


oTh .  grad  v 
a  a 


-  d(pkv  -  T.  pgv)  . 

a  aa 

6.  Constitutive  Equations  for  a  Special  Two-Component  Mixture 

Although  constitutive  equations  may  be  written  for  a  general  mixture  and  the 
restrictions  imposed  on  them  by  the  entropy  inequality  (5. 1 5)  may  be  found,  the  resulting 
algebraic  manipulations  are  quite  tedious.  I  consider  here  a  special 
two-component  mixture,  one  component  being  a  granular  material  as 
defined  by  Goodman  and  Cowin  (  2] ,  the  other  a  fluid  of  complexity  1.  In 
particular 


k  =  0 ,  k  =  0 . 


(6.1) 


. . . . . . . . . . 


Define 


L  -  {pcf  P  ,  grad  p  ,  x,  grad  x  .  a,  grad  a  ,  k,  v.  grad  v'.  v}  ,  (6.2; 

q  q  a  a  a  7Q  aa 


where  p  CK)  =  p fx  ,  O).  and  Q  =  I,  2  » 
cP  c  Q  Q 


Assume  that  each  cf  the  quantities 


4s  h>  2>  £>  1>  E1  *  cx,  c,  h,  v,  g 
a  o  a  a  a  o  a  0 


depends  on  L.  The  material  defined  by  this  constitutive  assumption  will 
be  called  a  mixture  of  a  granular  medium  with  a  fluid. 

Constitutive  equations  are  often  assumed  to  be  subject  to  certain 
restrictions,  one  of  which  is  called  the  "principle  of  fra  me -indifference". 
A  change  of  frame  is  defined  by 


jf  =  * "  +  9<z  '  *J.  QQT  =  l, 


p  -  p. 

£  S 


(6. 4; 


This  is  a  generalization  of  the  usual  definition. 


Define 


% 

Q  grad  $  ,  k,  v  ,  Q  grad  v,  v)  . 

2  o  o  a 


Consider  scalar-,  vector-,  and  tensor- valued  functions  s(L),  w(L), 

T(L)  .  These  functions  are  said  to  be  fra  me -indifferent  if 
o 

s(L* )  =  s(L)  , 

jy(L*)  =  Qw(L),  (6.5) 

TO* )  =  9T(L)9T  . 
o  o 

I  postulate  that  each  of  the  quantities  in  (6.3)  is  frame -indifferent.  Then 
by  a  familiar  argument,  dependence  of  these  functions  on  L  reduces  to 
dependence  on  L,  where 

L  =  {p  .  p  ,  grad  p,  v,  D,  S,  «,  grad  «,  k,  v,  grad  v,  v>,  (6.6) 

(TO  00  Z  a  00 

the  dependence  on  vectors  is  only  through  their  inner  products,  and 


a  =  W  -  W,  (6.7) 

i  2 

N  S 

v  =  x  -  x  =  u  -  u  .  (6.8) 

M  M  M  M  M 

12  12 
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Furthermore,  representation  theorems  for  (6.3)  as  functions  of  L  are  known, 
although  quite  complicated.  In  the  special  case  where  the  dependence  on 
the  vector  and  tensor  variables  in  L  is  linear, 

s(L)  =  s,  (6.9) 


w(L)  =  2 w  3raa  p  +  w^v  +  w,j,  grad  4  +  Ewv  grad  v  , 


(6.10) 


I(L)=-pl.  +  2  e  (tr  D)1  +  22  t,  D  , 
a  0  ab  b  ab  b 


(6.11) 


T(L)  =  -HO  , 
a 


(6.12) 


s  a 

where  T  and  T  are  the  symmetric  and  skew  parts  of  T ,  the  summations 
oo  _  o 

are  over  the  two  components,  and  the  functions  s,  w  ,  w  ,  w  ,  w  ,  p , 

op  D  1  aV  o 

i  r|  ,fi  are  functions  of  l ,  where 
ab  *  ob 


*  =  (P0>  P,  k>  v,  v) 

o  o  2  a  o 


(6.13) 


I  impose  the  restriction  that  L  and  I  do  not  in  fact  contain  p  , 

1° 

so  that  the  first  constituent  is  in  fact  a  fluid.  In  addition,  I  assume  that 


c  =  0, 
a 


(6.14) 


thus  excluding  certain  physical  phenomena,  including  chemical  reactions. 

By  (2.19),  (2. 22)3,  (5.19),  (6.7)  and  some  straightforward  algebra, 
‘  it  may  be  shown  that 


-10- 
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(6.15) 


2tr(T  G)  =  tr(T  D  +  TD  +  TfiT) 

ft  w  ^  W  «W  M  ' 

11  22  i 


o  o 


Also,  by  (2.23)  and  (6.8) 

w 


_  +  + 

2  m .  u  =  m  •  v  . 

•*  <•  ~  ~  > 

a  a  i 


(6.16) 


and  by  (2. 23)t 


4*  N  j.  x  ». 

2  v(v  -  v)  =  V(v  -  v)  . 

oo  112 


By  (1.6),  (1.10),  (2.11)  and  (6.8), 


(6.17) 


Epx*  u  = 
oa  a 


plp2 


v  •  v  . 


(6.18) 


It  is  easily  shown  that 


x  =  x  +  Gu . 

iv  «v  ^  7 

0  0  ft  o 


(6.19) 


so 


2  px  •  u  =  2  px  .  u  +  2  pii .  Gu  , 
oo  a  qq  a  aa  a  a 

< 

=  2"  pjc  •  u  +  2  pu*Du  , 

aa  o  aft  oo 


(6.20) 


or,  by  (6.18) 


2  px  •  u  = 

r  ««  «v 

0  0  ft  P 


plp2 


v  •  v  +  T  pu  •  Du  . 

^  r<V  M 

00  0  0 


(6.21) 
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By  a  set  of  steps  similar  to  those  leading  to  (6.  21),  it  may  be 
shown  that,  for  a  two- constituent  mixture  of  granular  materials 


„  "  V 

E(pkv  +  pkv)(v-  v) 
ooa  ooa  Q 

=  (v  -  v)  +  — — •  (kkv  -  kkv](v-v)  .  (6.22) 

Pk  1  2  Pk  121  122  1  2 


However,  in  the  special  case  k  =  0,  by  (1.11) 


V 


(6.23) 


and 


E  (pkv+  pkv)(v-v)  =  o  . 
qoq  aaa  a 

By  (6.14)  -  (6.17),  (6.21)  and  (6.24),  (5.28)  becomes 

$  •  ~  plp2  .  „  „  i 

pT1  ~  -  p^i>  div  $  -  *[-—  v  •  v  +  E  pu  .  Du] 

«  *  P  to  o  o 

N  +  +  '  # 

+  [g  -  ET  u  -  h(v-v)]  .  grad  $  +  $pm •  v  +  apv(v-v) 

ooll  1  11 

-  -  tr(T  - ET)(E pD  +  2u®  grad  p) 

P  0000  0 

s  s  a  t 

-  ■a  tr(TD  +  TD  +  Tfi  ) 

11  2  2  1 

-  $(h  -  Eh)  •  grad  v  -  4Eh  •  grad  v* 

o  2  0  o 


-  *{pkv2  -  r  pgv)  , 
so  a 


(6.24) 


(6.25) 


-12- 


4 


1  ?' 


wf\  fc-aiM  fcx  r 


A  result  derived  in  a  fashion  similar  to  (6 . 19 )  is 


4*  =  -2“^  (u  •  grad  p  +  p  tr  D) 

®P  a  a  a 

a 


"  Si(grad~p)  *  I  G  9rad  P  +  (grad  grad  P  }H 


+  (tr  D)grad  p  +  p  grad(tr  D)] 
a  a  a  a 

+  .  v  +  2  tr  b+  tr  Q  -4-  $ 

9v  ~  8D  ~  8fi  ~  8* 

w  *v  Q  ^ 

a 

+  -r . %~r:  •  17  grad  $  +  ^  k 

3(grad  fl)  dt  3k  _ 

2 

+  E  (v  -  u  •  grad  v) 


(6.31) 


+  S-— Vt  *  17  (grad  v 
3(grad  v)  dt 

a  11 


)  +  2“^-  •  ~(grad  v  ) 
dv  dt  *  Q 

a 


It  is  obvious  that  the  computation  of  div  g  will  involve  gradients 
of  second-order  tensors,  which  are  third-order  tensors.  Such  quantities 
are  in  general  ill -adapted  to  the  system  of  notation  used  here.  I  choose 
a  system  of  orthogonal  Cartesian  coordinates,  and  establish  the  notation 
used  in  terms  of  components  referred  to  these  coordinates. 

Let  a  be  a  vector  and  £  be  a  tensor.  Read  the  symbol  " 
as  "has  Cartesian  components",  so  that 


a  ~  a. 


<7  ~  a..  . 
~  jk 


Then  by  definition 


-H391 


>  i  ft}*  * ■  jkMLLila  ififerdUSlrfaliii 


9rad2  "V*  ’ 


where  the  comma  indicates  partial  differentiation,  and 


Define 


Note  the  identity 


da  da 


Jk 


trt^-arad 


fiiJ,k  =  (^ki,j  ‘  ^kj,i}  ’  ^ki,j  "  ^kj,i 


>, 


(6, 


which  will  be  abbreviated  as 


X  X 

grad  12  =  {grad  D)  ~  {grad  D )  . 

T  2 


{6. 


Differentiating  {6. 30), gives 


grad  v  s  —  grad  v+  {grad  grad  v)u 
0  0  {[  0 


{6. 


Computing  div  £  is  now  a  straightforward  task.  By  {6.35)  and  {6.36} 


34) 


35} 


.36) 


dlV  i  =  ^  ‘  grad  po  +  r  if  •  9r6  - 
2°  2  o  u 

+  Ttr  9{grad  p)  9rad  9rad  P  +  “  •  "  wrad  Y 

Q  ° 

+  2tr8D  9rad  5  +  tr^  (grad  Q)  -  tr ^  (grad  D* 


(6.37) 


grad  fl  +  tr 


i(5radl)  grad  9rad  * 


+  •  9rad  ’  ♦  *  *ii^r  «-  9rad  , 

2  -  °  «  0 

+  T  9v  *  fdt  grad  v  +  (grad  ^rad  v)u]  . 


Note  that,  by  (2.12)  (6.1)  and  (6.14) 


pk  s  -div(pku), 
222 


so  that 


pk  =  -k  dlv(pu)  «  pu  •  cirad  k 
2  tz  ZZ  ’  2 


By  (6.23)  and  (6. 391 


pkv2  =  -kv2  div(pu)  -  pv2u  •  grad  k  . 
22  22  22  2  2 


Also 


k  =  -  u  •  grad  k  . 
2  2  2 


Substituting  (6.33),  (6,40)  and  (6.41)  into  (6.25)  then  yields 


(6.38) 


(6.39) 


(6.40) 


(6.41) 
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A  thermodynamic  process  is  defined  in  a  fashion  similar  to  that  for  a  conven¬ 
tional  continuum.  Each  balance  law  with  the  exception  of  the  one  for  mass  contains 
one  term  which  may  be  adjusted  arbitrarily,  and  there  exists  at  least  one 

thermodynamic  process  in  which  v  ,  D ,  grad  D,  Q  ,  grad  p  , 
d  o  ,  s  Z  a° 

grad  grad  p  ,  --j^grad  grad  9rad  ^  v,  --grad  v  ,  g*ad  grad  v 

o  a  a  q 

and  grad  k  may  be  chosen  arbitrarily  and  independently  of  any 
2 

jjk 

other  term  in  the  inequality.  This  implies  the  following  results  : 


a.  4>  is  independent  of  D,  n,  grad  4,  v  f 

a  (x 


b. 


<Hi  P1P2 
8v  2  ~  * 


so  that 


(6,43) 


-  * 

*  =  »Mp0,  P,  9rad  p,  i,  k,  V,  grad  v)  +  -±£  -  u-  u  .  (6.44) 

2  0  a  20  a^pQtt 

Furthermore,  4*j  depends  on  grad  p  and  grad  v  only  through  their  inner 

a  a 

products  grad  p  •  grad  p  ,  grad  p  •  grad  v  ,  grad  v  »  grad  v  . 

a  b  a  b  0  b 


c. 


isk  _  JL 

a*  ~  .2 


ak 


d.  =  P«uH£  -  v  ")  . 


/M 

'ak 


e.  (• 


ii. 


8(grad  a) 


)S  =2  , 


(6.45) 

(6.46) 


(6.47) 


These  results  reduce  to  those  of  Muller  [  3]  in  the  case  where  both 
continua  are  ordinary  fluids.  I  have  corrected  a  minor  printer’s  error  in 
result  b. 
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5 

where,  as  before,  the  superscript  denotes  the  symmetric  part  of 
the  indicated  tensor. 


f-  =2-  (6.48)- 

0 

Here  (6.44)  has  been  used. 

g.  Again,  I  revert  to  Cartesian  components. 


(pp* 


3vji 


I 


1  9P,k 

1 


U 


94>l 

“Vi 

4, 


f+V 

an 


)i 


o 

8?ii  ’ 


(6.49) 


(pp2*  6«+«Vi,"«v'8|; 

s 


8iJ 

)  3  =  o. 


(6.50) 


where  the  superscript  denotes  the  symmetric  part  with  respect  to 


i  and  J. 


d$  9^I 

h.  «h  =  -^  +  flp  Q/  *  , 

7  3V  r  3 (grad  v)  * 

1  1  1 

(6.51) 

5  W  P  8(grad  v,  l2z- 

(6.52) 

l-  Wadv,  (4  -P'+iHl* 

(6.53) 

a 


where  in  deriving  (6.  53),  (6.  51)  and  (6.  52)  have  been  used. 

j.  *  is  independent  of  grad  p  . 

*  _o 
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raspa 


I  note  that 


Furthermore 


grad  v  =  D  -  D  +  fit  . 

1  2 


div  pu  =  p  div(x  -  x)  +  u  •  grad  p  , 

22  2  2  ~  2  2 

=  ptrD-ptrD+u*  grad  p  . 
2  2  2  ~  2  2 

However,  by  (5.22) 


p  tr  D  =  p  tr  D  +  p  tr  D  +  u  '  grad  p.  f  u  *  grad  p  . 

r  z  2  r  1  i  2 

Substituting  *  .56)  into  (6.55)  and  taking  (1.6)  into  account,  I 
obtain 

pip2  S 

i  v  p  u  = -  (tr  D  -  tr  D)  -  —  u  *  grad  p 

22  p  2  1  PT  1 

P 

+  -  u  •  grad  p  . 

p  :>  2 

Substitution  of  (6.54),  and' (6. 57)  into  (6.42)  yields  the 
i  .sidual  inequality 


(6.54) 


(6.55) 


(6.56) 


(6.57) 
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wmm  ^pill'^WFPPBIW'lllfllHiiPJi'H.11:  \  m  IWV  MW*  il 111 i  up 


The  consequences  of  resuits  a.  -  f.  may  be  investigated  by 

standard  methods  and  constitute  generalizations  of  known  results  {see, 

e.g. ,  [  3])  for  a  mixture  of  two  constituents  to  the  case  where  one  of 

the  components  is  a  granular  material.  In  particular  a.  and  b.  indicate 

that  the  Helmholtz  free  energy,  in  addition  to  being  independent  of  peculiar 

stretchings  D,  the  relative  spin  g,  and  the  coldness  gradient,  grad  a,  is 
a 

independent  of  the  peculiar  volume  distribution  velocities  v  .  Furthermore 

a 

its  dependence  on  diffusion  velocities  is  made  explicit  by  (6.44).  Result 
c.  is  a  familiar  relation  from  thermodynamics,  but  result  d.  is  new. 
Results  h.  generalize  analogous  results  of  Goodman  and  Cowin  indicating, 
as  might  be  expected,  a  contribution  tc  peculiar  equilibrated  stress  due  to 
diffusion.  Result  g.  constitutes  a  set  of  restrictions  on  the  forms  of 
the  entropy  flux  and  Helmholtz  free  energy.  The  restrictions  imposed 
by  the  familiar  relations  e.  and  f. ,  and  the  new  relation  i.  are  easily  made 
explicit  by  a  method  due  tc  Muller.  Solving  these  equations,  I  obtain,  in 
Cartesian  tensor  notation, 


♦i  ■  *fiVi  =  ,(*>i  pri  > + >  prj  wlfj f 

+  tijkpj  rk  +  tij  +  ’ 


16.59) 
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(6.60) 


*1  +  =  (^ijkl  \,i  *,k  +  ^ijk  P,J  +I^ijk  *,J  + 


+  >  +  ^  V,J  +  r9iJP7,J  +  Aoi* 


Here  the  coefficients  A,.  r,  A:\are  skevt*~in.ail  indices,  are  independent 


of  grad  and  a  subscript  1  or  2  indicates  independence  of  grad  p 

1 

and  grad  v  or  grad  p  and  grad*,  respectively. 

12  2 
It  follows  that 


P&PrVi  = 

1  2 


«V2>k 


+  X-1+r21!2-)  +  A2l) 


+  A 
1 


«v1,‘+ri«p1,)  +  ^) 


1 


(6.61) 


Equations  (6.59)  -  (6.61)  may  be  used  to  derive  a  more  explicit 

form  for  £  by  straightforward  addition  and  substitution.  Furthermore, 

the  functions  A,  A,  r  may  be  recognized  as  combinations  of  the  derivatives 

of  pfli^v  with  respect  to  grad  p,  grad  v  and  the  values  of  these 

a  q 

-  derivatives  with  one  or  more  of  the  parameters  vanishing.  Finally,  it 

A 

may  be  recalled  that  <|iT  and  §  depend  on  grad  p  ,  grad  v  and  grad  $ 

1  o  a 

only  through  their  inner  products,  thus  yielding  an  alternative  form  for  ipj. 

All  of  these  results  involve  straightforward  calculations  and  yield  somewhat 
complicated  results.  Since  they  are  not  needed  in  their  full  generality  in 
the  following  sections,  they  are  not  recorded  here. 
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7 .  Equilibrium 


X  =  {grad  p,  D ,  g  ,  grad  4 ,  v  ,  grad  v  ,  v}  . 


(7.1) 


a  a 


o  a 


If  X.  =  0,  the  mixture  is  said  to  be  in  equilibrium.  Let  the  entropy 
production  cr  be  the  left-hand  side  of  (6.58),  so  that  inequality  is 


<r  >  0  . 


(7.2) 


Then  o-  has  a  minimum  at  equilibrium.  Necessary  condtions  for  this  are 


7x7  =  °  at  V°* 

A 


(7.3) 


8X  ax.  I  is  non -negative  definite  at  X^  =  0  .  (7.4) 

A  B  i 

By  (2.25)^  and  (6.59)  -  (6.61),  consequences  of  (7.3)  are 


!*  =  o, 

dp  * 


(7.5) 


AT  =  App  1  -  ^  m  -  4p  1  , 
~  8p  ~  9v  ~  ’ 


1  1 


(7.6) 


AT  =  App  1  +  ^  s  -flpl  , 

2  2  9p  ~  9~  2 


a 

T  =  0  , 
1 


(7.7) 


(7.8) 


3  =  2 


(7.9) 
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m  =  0  , 

(7.10) 

1 

4 

i5  =°> 

(7.11) 

a 

pg  = 
11 

8dj  + 

-ps? " > 

i  i 

(7.12) 

pg  - 
22 

‘p  8?  + 

2 

(7.13) 

where  all  quantities  are  evaluated  at  equilibrium. 

Equations  (7.5)  -  (7.10)  are  familiar  results  from  the  thermostatics 
of  a  mixture  of  two  ordinary  fluids,  and  have  been  commented  upon 
extensively  by  Muller  [  3] .  Unlike  the  single  component  granular  material 
of  Goodman  and  Cowin,  a  mixture  of  a  granular  material  and  a  fluid  is 
incapable  of  sustaining  a  shear  stress  in  equilibrium.  The  results  (7.12) 
and  (7.13)  are  analogous  to  those  obtained  in  elasticity  theory,  stating 
that  the  intrinsic  body  forces  are  derived  from  the  Helmholtz  free  energy 
and  the  growth  of  equilibrated  force. 

8.  A  Linear  Theory 

Consider  an  expansion  about  equilibrium,  linear  in  X.  as  defined 

7  A 

by  (7.1)  of  the  quantities  4>,  <£  ,  T,  m  , ,  h  ,  v  ,  g  . 

o  1  o  1  o 

By  results  a.  and  b.  in  Section  6 

=  4»(P0>  P,  grad  p  ,  «,  k,  v  ,  grad  v )  .  (8.1) 

2Q  o  2  a  a 
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However,  by  the  assumption  of  frame- indifference,  ^  depends  on  grad  p 

>) 

and  grad  v  through  thsir  inner  products,  which  are  nonlinear.  Thus 


=  t>(pn,  9,  P>  k,  v>  v)  • 

2  12  2  12 


(8.2) 


By  (6.10),  linear  representations  for  q,  m,  and  are 

J 

%  =  "Kp  grad  *  2'k  grad  p-  k  v  -  X  grad  v  f  (8.3) 

*  (jP  a  ^  $  a 

-f 

m  =  -m  grad  3  -  Xm  grad  p  -  m_v  -  x  m  grad  v  ,  (8.4) 

T  T  a?  o  ^  av  a 

$  =  -1^  grad  a  -  X  K  grad  p  -  K^v  -  SKy  grad  V ,  (8.4) 


v/here  the  coefficients  are  functions  of  p  ,  p,  k,  »>,  v.  The  represents  - 

2°  o  2oo 

tion  (8.5)  for  4>  is  further  restricted  by  (6.47),  (6.48),  (6.53)  and  (8.2). 
Straightforward  computation  leads  to  the  conclusion  that  K_,  K  ,  K 

1  D  V 
0K  0 


vanish,  giving 


W 


(8.6) 


In  this  case  (6.49)  and  (6.50)  are  satisfied  identically.  It  is  seen  that, 
unlike  the  special  case  of  a  single  granular  material  of  Goodman  ano 
Cowin  [  2],  the  entropy  flux  does  not  have  the  classical  form  of  the  heat 
flux  multiplied  by  the  coldness.  It  is  in  addition  affected  by  the  diffusion 
velocity,  with  coefficient  depending  on  density,  volume  distribution  and 
its  velocity,  coldness,  and  equilibrated  inertia. 
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v-4" 


l 


1 


S'4" 


K  >0 
T  -  ' 


«  =  0,  K  =0, 

2P 


1* 


6Kt 


_ p  dip 

$om  -  -  — —  -  3p  r*  , 

•  ;p  dp  ,  3p 

i  2  i 

3pm  =  -  +  3p 

,P  3p  j  i>P 

w  *>  *  ws 


(8.13) 


t  >  0 

V_  u} 

i  >o, 

1 

1 

11 

22 

11 

22 

h  >0, 

h  >of 

T| 

»} 

11 

22 

11 

22 

J. 

4 


12  Zi 


12  21 


K-  >0, 


2^mD  >  0  , 


8K 


P  P 


3K. 


n  *»  i  v**p'  o  i 

2«T(i8pmD)  >  (—  +  «D  4  *pmT  p  -j  -K-j^+  “D  +  ifiij  +  7-) 


P  P 

i 

P 


Finally,  the  determinant 
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n 

4  +  4 

X  +  6  -  pp 

X  4  6  4  pp 

u 

12  21 

11  11  1 

12  21  1 

4  +  4 

24 

X  4  6  -  pp 

X  4  64  pp 

12  21 

22 

21  12  2 

22  22  2 

X  +  6  -  pp 

X  +  6  -  pp 

2(4  -  pu) 

4  A  4  +  p(u-u) 

11  11  1 

21  12  2 

11  1 

12  21  2  1 

X  4  6  4  pp 

X  4  5  4  pp 

4  +  4  +  p(«  -  t> 

2(4  +  pu) 

12  1 

22  22  2 

12  21  21 

22  2 

(8.14) 


must  be  positive  semi-definite.  The  restrictions  (8.13)  are  classical 
restrictions  for  a  mixture  of  two  fluids  as  given  by  Muller.  The 
inequalities  (8.14)  include  (8.13)?  but  are  otherwise  new  to  this 
theory. 


The  form  of  the  entropy  flux  <{>  ir.  the  linear  theory  can  be  found. 
By  (8.13), 


%2 


P  P 
2 


KD  =  )  +  (7  p  -  7  P  )  ]  . 


12 


dp  dp 


(8.15) 


1  ^  2 


1  2 

Substituting  (6.  23),  (7,6)  and  (8.6)  into  (5.12),  and  comparing  with  (8.15) 
substituted  into  (8.6)  gives 


8K 

4>  =  flg  -  3v[pp  ("^  -  ^)  -  — ^  (v  -  v)l  . 

~  »  ~  l?  3P  9P  dv  \  ,J 

1  2  l  C 


(8.16) 


This  generalizes  Muller's  (7.21). 


#1391 


-29- 


ACKNOWLEDGEMENT 


My  thanks  go  to  the  Mathematics  Research  Center  and  the  Soil 
Scienca  Department  of  the  University  of  Wisconsin  for  support  of 
this  work. 


-30- 


#1391 


REFERENCES 


1.  S.  L.  Passman:  On  the  balance  equations  for  a  mixture  of 
granular  materials.  Mathematics  Research  Center,  University 
of  Wisconsin-Madison,  TSR  #1390  (1973). 

2.  M.  A.  Goodman  and  S.  C.  Cowin:  A  continuum  theory  for  granular 
materials .  Archive  for  Rational  Mechanics  and  Analysis  44  , 

4,  249-266  (197  2). 

3.  I.  Muller:  A  thermodynamic  theory  of  mixtures  of  fluids. 

Archive  for  Rational  Mechanics  and  Analysis  28,  1-39  (1968). 


I 


K 


#1391 


-31- 


